webGQT: A Shiny Server for Genotype Query Tools for Model-Based Variant Filtering

Genotype Query Tools (GQT) were developed to discover disease-causing variations from billions of genotypes and millions of genomes, processes data at substantially higher speed over other existing methods. While GQT has been available to a wide audience as command-line software, the difficulty of constructing queries among non-IT or non-bioinformatics researchers has limited its applicability. To overcome this limitation, we developed webGQT, an easy-to-use tool with a graphical user interface. With pre-built queries across three modules, webGQT allows for pedigree analysis, case-control studies, and population frequency studies. As a package, webGQT allows researchers with less or no applied bioinformatics/IT experience to mine potential disease-causing variants from billions. Results webGQT offers a flexible and easy-to-use interface for model-based candidate variant filtering for Mendelian diseases from thousands to millions of genomes at a reduced computation time. Additionally, webGQT provides adjustable parameters to reduce false positives and rescue missing genotypes across all modules. Using a case study, we demonstrate the applicability of webGQT to query non-human genomes. In addition, we demonstrate the scalability of webGQT on large data sets by implementing complex population-specific queries on the 1000 Genomes Project Phase 3 data set, which includes 8.4 billion variants from 2504 individuals across 26 different populations. Furthermore, webGQT supports filtering single-nucleotide variants, short insertions/deletions, copy number or any other variant genotypes supported by the VCF specification. Our results show that webGQT can be used as an online web service, or deployed on personal computers or local servers within research groups. Availability webGQT is made available to the users in three forms: 1) as a webserver available at https://vm1138.kaj.pouta.csc.fi/webgqt/, 2) as an R package to install on personal computers, and 3) as part of the same R package to configure on the user's own servers. The application is available for installation at https://github.com/arumds/webgqt.Peer reviewe

Binary Interval Search (BITS): A Scalable Algorithm for Counting Interval Intersections

Motivation: The comparison of diverse genomic datasets is fundamental to understanding genome biology. Researchers must explore many large datasets of genome intervals (e.g., genes, sequence alignments) to place their experimental results in a broader context and to make new discoveries. Relationships between genomic datasets are typically measured by identifying intervals that intersect: that is, they overlap and thus share a common genome interval. Given the continued advances in DNA sequencing technologies, efficient methods for measuring statistically significant relationships between many sets of genomic features is crucial for future discovery. Results: We introduce the Binary Interval Search (BITS) algorithm, a novel and scalable approach to interval set intersection. We demonstrate that BITS outperforms existing methods at counting interval intersections. Moreover, we show that BITS is intrinsically suited to parallel computing architectures such as Graphics Processing Units (GPUs) by illustrating its utility for efficient Monte-Carlo simulations measuring the significance of relationships between sets of genomic intervals

The structural variation landscape in 492 Atlantic salmon genomes

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species.Peer reviewe

Combating subclonal evolution of resistant cancer phenotypes

Metastatic breast cancer remains challenging to treat, and most patients ultimately progress on therapy. This acquired drug resistance is largely due to drug-refractory sub-populations (subclones) within heterogeneous tumors. Here, we track the genetic and phenotypic subclonal evolution of four breast cancers through years of treatment to better understand how breast cancers become drug-resistant. Recurrently appearing post-chemotherapy mutations are rare. However, bulk and single-cell RNA sequencing reveal acquisition of malignant phenotypes after treatment, including enhanced mesenchymal and growth factor signaling, which may promote drug resistance, and decreased antigen presentation and TNF-α signaling, which may enable immune system avoidance. Some of these phenotypes pre-exist in pre-treatment subclones that become dominant after chemotherapy, indicating selection for resistance phenotypes. Post-chemotherapy cancer cells are effectively treated with drugs targeting acquired phenotypes. These findings highlight cancer's ability to evolve phenotypically and suggest a phenotype-targeted treatment strategy that adapts to cancer as it evolves

Mapping and characterization of structural variation in 17,795 human genomes.

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing

A Soluble Acetylcholinesterase Provides Chemical Defense against Xenobiotics in the Pinewood Nematode

The pinewood nematode genome encodes at least three distinct acetylcholinesterases (AChEs). To understand physiological roles of the three pinewood nematode AChEs (BxACE-1, BxACE-2, and BxACE-3), BxACE-3 in particular, their tissue distribution and inhibition profiles were investigated. Immunohistochemistry revealed that BxACE-1 and BxACE-2 were distributed in neuronal tissues. In contrast, BxACE-3 was detected from some specific tissues and extracted without the aid of detergent, suggesting its soluble nature unlike BxACE-1 and BxACE-2. When present together, BxAChE3 significantly reduced the inhibition of BxACE-1 and BxACE-2 by cholinesterase inhibitors. Knockdown of BxACE-3 by RNA interference significantly increased the toxicity of three nematicidal compounds, supporting the protective role of BxACE-3 against chemicals. In summary, BxACE-3 appears to have a non-neuronal function of chemical defense whereas both BxACE-1 and BxACE-2 have classical neuronal function of synaptic transmission

Human-based fibrillar nanocomposite hydrogels as bioinstructive matrices to tune stem cell behavior

The extracellular matrix (ECM)-biomimetic fibrillar structure of platelet lysate (PL) gels along with its enriched milieu of biomolecules has drawn significant interest in regenerative medicine applications. However, PL-based gels have poor structural stability which severely limits its performance as a bioinstructive biomaterial. Here, rod-shaped cellulose nanocrystals (CNC) are used as a novel approach to modulate the physical and biochemical microenvironment of PL gels enabling their effective use as injectable human-based cell scaffolds with a level of biomimicry that is difficult to recreate with synthetic biomaterials. The incorporation of CNC (0 to 0.61 wt.%) into the PL fibrillar network during the coagulation cascade leads to decreased fiber branching, increased interfiber porosity (from 66 to 83%) and modulate fiber (from 1.4 ± 0.7 to 27 ± 12 kPa) and bulk hydrogel (from 18 ± 4 to 1256 ± 82 Pa) mechanical properties. As result of these physicochemical alterations, nanocomposite PL hydrogels resist to the typical extensive clot retraction (from 76 ± 1 to 24 ± 3 at Day 7) and show favored retention of PL bioactive molecules. The feedback of these cues on the fate of human adipose-derived stem cells is evaluated, showing how it can be explored to modulate the commitment of encapsulated stem cells toward different genetic phenotypes without the need for additional external biological stimuli. These fibrillar nanocomposite hydrogels allow therefore to explore the outstanding biological properties of human-based PL as an efficient engineered ECM which can be tailored to trigger specific regenerative pathways in minimal invasive strategies.The authors thank the Hospital da Prelada (Porto, Portugal) for providing adipose tissue samples. The authors acknowledge the financial support from project Recognize (UTAP-ICDT/CTM-BIO/0023/2014), project NORTE-01-0145FEDER-000021 supported by the Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF), the European Union Framework Programme for Research and Innovation HORIZON 2020, under the TEAMING Grant agreement No. 739572 – The Discoveries CTR EU, Forecast 668983, Marie Skłodowska-Curie grant agreement No. 706996 (PrinTendon) and CHEM2NATURE 692333; FCT/MCTES (Fundação para a Ciência e a Tecnologia/ Ministério da Ciência, Tecnologia, e Ensino Superior) and the Fundo Social Europeu através do Programa Operacional do Capital Humano (FSE/POCH) in the framework of PhD grant PD/59/2013 – PD/BD/113807/2015 for BBM, Post-Doc grant SFRH/BPD/112459/2015 for R.D.info:eu-repo/semantics/publishedVersio

Nuclear Scaffold Attachment Sites within ENCODE Regions Associate with Actively Transcribed Genes

The human genome must be packaged and organized in a functional manner for the regulation of DNA replication and transcription. The nuclear scaffold/matrix, consisting of structural and functional nuclear proteins, remains after extraction of nuclei and anchors loops of DNA. In the search for cis-elements functioning as chromatin domain boundaries, we identified 453 nuclear scaffold attachment sites purified by lithium-3,5-iodosalicylate extraction of HeLa nuclei across 30 Mb of the human genome studied by the ENCODE pilot project. The scaffold attachment sites mapped predominately near expressed genes and localized near transcription start sites and the ends of genes but not to boundary elements. In addition, these regions were enriched for RNA polymerase II and transcription factor binding sites and were located in early replicating regions of the genome. We believe these sites correspond to genome-interactions mediated by transcription factors and transcriptional machinery immobilized on a nuclear substructure

Impact of COVID-19 on cardiovascular testing in the United States versus the rest of the world

Objectives: This study sought to quantify and compare the decline in volumes of cardiovascular procedures between the United States and non-US institutions during the early phase of the coronavirus disease-2019 (COVID-19) pandemic. Background: The COVID-19 pandemic has disrupted the care of many non-COVID-19 illnesses. Reductions in diagnostic cardiovascular testing around the world have led to concerns over the implications of reduced testing for cardiovascular disease (CVD) morbidity and mortality. Methods: Data were submitted to the INCAPS-COVID (International Atomic Energy Agency Non-Invasive Cardiology Protocols Study of COVID-19), a multinational registry comprising 909 institutions in 108 countries (including 155 facilities in 40 U.S. states), assessing the impact of the COVID-19 pandemic on volumes of diagnostic cardiovascular procedures. Data were obtained for April 2020 and compared with volumes of baseline procedures from March 2019. We compared laboratory characteristics, practices, and procedure volumes between U.S. and non-U.S. facilities and between U.S. geographic regions and identified factors associated with volume reduction in the United States. Results: Reductions in the volumes of procedures in the United States were similar to those in non-U.S. facilities (68% vs. 63%, respectively; p = 0.237), although U.S. facilities reported greater reductions in invasive coronary angiography (69% vs. 53%, respectively; p < 0.001). Significantly more U.S. facilities reported increased use of telehealth and patient screening measures than non-U.S. facilities, such as temperature checks, symptom screenings, and COVID-19 testing. Reductions in volumes of procedures differed between U.S. regions, with larger declines observed in the Northeast (76%) and Midwest (74%) than in the South (62%) and West (44%). Prevalence of COVID-19, staff redeployments, outpatient centers, and urban centers were associated with greater reductions in volume in U.S. facilities in a multivariable analysis. Conclusions: We observed marked reductions in U.S. cardiovascular testing in the early phase of the pandemic and significant variability between U.S. regions. The association between reductions of volumes and COVID-19 prevalence in the United States highlighted the need for proactive efforts to maintain access to cardiovascular testing in areas most affected by outbreaks of COVID-19 infection